Observation of an Inverse Energy Cascade 
in Developed Acoustic Turbulence in Superfluid Helium 
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We report observation of an inverse energy cascade in second sound acoustic turbulence in Hell. 
Its onset occurs above a critical driving energy and it is accompanied by giant waves that constitute 
an acoustic analogue of the rogue waves that occasionally appear on the surface of the ocean. The 
theory of the phenomenon is developed and shown to be in good agreement with the experiments. 
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A highly excited state of a system with numerous de- 
grees of freedom, characterized by a directional energy 
flux through frequency scales, is referred to as turbulent 
P, 0| • Like the familiar manifestations of vortex turbu- 
lence in fluids, turbulence can also occur in systems of 
waves, e.g. turbulence of sound waves in oceanic waveg- 
uides Q, magnetic turbulence in interstellar gases [jj, 
shock waves in the solar wind and their coupling with 
Earth's magnetosphere [5j, and phonon turbulence in 
solids Q. Following the ideas of Kolmogorov, the uni- 
versally accepted picture says that nonlinear wave inter- 
actions give rise to a cascade of wave energy towards 
shorter and shorter wavelengths until, eventually, it be- 
comes possible for viscosity to dissipate the energy as 
heat. Experiments and calculations show that, most of 
the time, the Kolmogorov picture is correct 0, @, @] • 

We demonstrate below that this picture is incom- 
plete. Our experiments with second sound (temperature- 
entropy) waves in Hell show that, contrary to the con- 
ventional wisdom, acoustic wave energy can sometimes 
flow in the opposite direction too. We note that inverse 
energy cascades are known in 2-dimensional incompress- 
ible liquids and Bose gases , and have been considered 
for quantized vortices [loj ]. 

We find that energy backflow in our acoustic system 
is attributable to a decay instability (cf. the kinetic in- 
stability in turbulent systems [ill]), controlled mainly by 
nonlinear decay of the wave into two waves of lower fre- 
quency governed by the energy (frequency) conservation 
law [2j u>i = ^2 + <-^3- Here uii — u-i^ki is the frequency 
of a linear wave of wave vector ki and U20 is the second 
sound velocity at negligibly small amplitude. The in- 
stability manifests itself in the generation of subharmon- 
ics. A quite similar parametric process, due to 4-wave 
scattering (modulation instability), is thought to be re- 
sponsible for the generation of large wind-driven ocean 



waves [121 ] . Decay instabilities (especially threshold and 
near-threshold behaviour) have been studied for e.g. spin 
waves [13(1 , magnetohydrodynamic waves in plasma 14(, 
and interacting first and second sound waves in superfluid 
helium near the superfluid transition [15j . 

We now discuss what happens to a system of acoustic 
waves far beyond the decay threshold. Modelling the re- 
sultant nonlinear wave transformations in the laboratory 
is a potentially fruitful approach that has already yielded 
important results for e.g. the turbulent decay of capillary 
waves on the surface of liquid H2 Here, we exploit 

the special properties of second sound [13] , which enable 
fundamental wave processes to be studied under labora- 
tory conditions. Its velocity u 2 depends strongly on its 
amplitude ST and can be approximated as 



u 2 = ?i2o(l + a5T), 



(1) 
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where the nonlinearity coefficient a [18j, which deter- 
mines the strength of the wave interactions, can be made 
large if the temperature is set right. 

Our experiments make use of the high Q cryoacousti- 
cal resonator described previously [7( , excited close to one 
of its resonant frequencies. It enables very large second 
sound standing wave amplitudes to be attained for mod- 
est levels of excitation. Its Q-factor, determined from the 
widths of longitudinal resonances at small heat flux den- 
sities, was Q ~ 1000 for resonant numbers p < 10 and 
Q ~ 3000 for 30 < p < 100. Fig.[T]presents typical results 
obtained when driving at a relatively high resonant fre- 
quency u>d (the 96th longitudinal resonance of the cell). 
Those in Figs. []Ja),(b) reproduce our earlier observation 
of the direct Kolmogorov-like cascade of second sound 
waves in He II Q, when driving on resonance. Figs.QJc)- 
(d) show that tiny changes in driving frequency can pro- 
duce marked changes in the shape and spectrum of the 
standing wave. The formation of the spectral peaks near 
i , i and I x uJd satisfies the frequency conservation law 
with u>i — ujd, supporting our inference that the insta- 
bility is controlled mainly by 3-wave interactions. When 
the instability develops, huge distortions of the initially 
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FIG. 1: (Color online) Evolution of the observed wave shape in 
the resonator (left column) and of the power spectrum of second 
sound standing waves (right) with increasing drive frequency Ud 
near the 96th resonance: ui^/2tt = 9530.8 Hz (a),(b); and 9535.2 
Hz (c),(d). The AC heat flux density was W = 42 mW/cm 2 . The 
temperature T = 2.08 K corresponded to negative nonlinearity. 
The fundamental and first harmonic in (b),(d) are indicated by 
vertical (green and blue) arrows; the low-frequency domain where 
the subharmonics appear are indicated by the horizontal range in 
(d) with vertical (red) arrows. The horizontal arrows in (a),(c) 
indicate the fundamental period of a wave at the driving frequency. 



periodic signal occur. Although it remains nearly peri- 
odic, its characteristic period exceeds that of the driving 
force and its amplitude can become more than twice that 
at the driving frequency. 

To characterise the instability quantitatively, we use 
the energy contained in the low-frequency part of the 
spectrum u> < u>d, 



Elf — 



dC_ 
df 



UJ<U) d 



(2) 



as an indicator. Fig. [2] shows the dependence of Elf 
on the AC heat flux density W, when driving on the 
96th resonance close to 2.08 K. For small W we did not 
observe any subharmonic generation at all 0]. Then, 
above a critical flux W c , E L f rose rapidly, suggesting 
that the phenomenon is of a threshold character. At 
T = T a = 1.88 K for which a vanishes [IH, no subhar- 
monics were observed, regardless of the magnitude of W, 
thus confirming the crucial importance of nonlinearity. 
For W above 10.4 mW/cm 2 , we observed a distortion 
of the signal similar to that shown in Fig. [He) and the 
formation of a few subharmonics. Further increase of 
W above 20 mW/cm 2 led to the generation of multiple 
subharmonics. These phenomena appear in the regime 
where the energy cascade towards the high frequency do- 
main (i.e. direct cascade, with a Kolmogorov-like spec- 
trum is already well-developed; see also Fig. Ill 

All the above results correspond to steady-state 
regimes, of the wave system. In Fig. [3] we illustrate 



12 

E 

8 



x 

Ll 

LU" 



1 1 


s - 


- -• 

-' \ 




✓ 




/ 

/ 

W '• 

1 ' 


"bio 


/v. 

A 

A A 
A 


- b 






/ 

/ 




9540 9544 


/ 




co d /27i (Hz) 



20 



40 

W (mW/cm 2 ) 



60 



FIG. 2: (Color online) The energy Elf contained in the low- 
frequency part of the spectrum as a function of the AC heat flux 
density W, while driving near to the 96th resonance for T ~ 2.08 K. 
The threshold value of W, marked by the (green) arrow, was W c = 
10.4 mW/cm 2 . The points are from experiment; dashed lines are 
guides to the eye. Inset: the dependence of E^p on uj^, measured 
for W = 55.6 mW/cm 2 ; the (red) arrow labels the maximum value 
of Elp , which is taken to the main figure. 



the transient processes observed after a step-like shift of 
the driving frequency from a frequency initially set far 
from any resonance to the 96th resonance frequency for 
W = 42.1mW/cm 2 , T = 2.08 K. We find that harmon- 
ics of the drive in the high-frequency spectral domain are 
formed very quickly, but that formation of the subhar- 
monics takes much longer: it took ~ 0.5 s here, and can 
reach several tens of seconds under some conditions [li| . 
It is evident from the inset in Fig. [3] that, as the in- 
stability develops, isolated "rogue waves" appear in the 
signal. As time evolves, the rogue waves appear more 
frequently and, at the later stages, they merge result- 
ing in the strong low-frequency modulation of the signal 
observed in the steady-state measurements (FigQJ. 

When the subharmonics appear, a marked reduction 
occurs in the energy contained in the high-frequency 
spectral domain: see Fig. 01 A substantial proportion 
of the energy then flows from the driving frequency u)d 
towards the low frequency domain lu < lu^ leading to 
an accumulation of wave energy there and a correspond- 
ing increase in wave amplitude. The reduction of wave 
amplitude seen in the high frequency spectral domain is 
indicative of the onset of energy backflow towards lower 
frequencies, i.e. a sharing of the flux between the direct 
and inverse energy cascades. The decrease in energy 
at high frequencies energy in 0.397s < t < 1.3 s is at- 
tributable to relaxation processes in the direct cascade. 
Redistribution of wave energy due to sharing of the en- 
ergy flux between the direct and inverse cascades starts 
at t = 1.3 s. Note that the transient evolution shown 
in Figs 3 and 4 is incomplete: the forward and inverse 
energy fluxes are still changing at t = 2s, leading us to 
anticipate further relaxation oscillations at longer times 
(inaccessible with our present equipment): the transient 
dynamics is clearly complex, and a full characterisation 



3 




FIG. 3: (Color online) Transient evolution of the 2nd sound 
wave amplitude ST after a step-like shift of the driving fre- 
quency to the 96th resonance at time t = 0.397 s. Signals in 
frames 1 and 3 are similar to those obtained in steady-state 
measurements, Fig. 1 (a) and (e) respectively. Formation of 
isolated "rogue" waves are clearly evident. Inset: example of 
a rogue wave, enlarged from frame 2. 
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FIG. 4: (Color online) Instantaneous spectra in frames 1 and 
3 of Fig. [3] The lower (blue) spectrum, for frame 1, shows the 
direct cascade only; the upper (orange) spectrum, for frame 3, 
shows both the direct and inverse cascades. The (green) arrow 
indicates the the fundamental peak at the driving frequency. 
Inset: evolution of the wave energy in the low- frequency and 
high frequency domains is shown by the (orange) squares and 
(blue) triangles respectively; (black) arrows mark the posi- 
tions of frames 1 and 3. 



will require further work. Absorption of wave energy at 
low frequencies is probably attributable to viscous drag 
of the normal fluid component on the resonator walls, 
given that bulk second sound damping is negligibly small 
in this frequency range: this would be consistent with the 
observed strong decrease of the resonator Q-factor below 
3 kHz. We observe hysteresis between increasing and de- 
creasing frequency scans (bars on data points in Fig. [5]). 
Its width, i.e. the region where the low frequency sound 
waves are metastable, was less than the viscous width of 
the resonance. 

To seek a more detailed understanding of wave energy 
transformation in acoustic turbulence, we used a tech- 
nique [19l | for direct numerical integration of the 2-fluid 
thermohydrodynamical equations [17], expanded up to 
quadratic terms in the wave amplitude. It representats 



the second sound waves in terms of Hamiltonian vari- 
ables @, as in earlier studies fgj of acoustic turbu- 
lence. Wave damping was, however, taken explicitly into 
account at all frequencies, a feature that is of key impor- 
tance for a correct description of subharmonic generation. 
The main results are summarised in Fig. 

It is evident (inset of Fig. [5]) that, for sufficiently high 
driving amplitude W, the wave develops an instability 
with respect to generation of low- frequency subharmonics 
of the driving force at uid- For zero detuning from a cavity 
resonance, the onset of the instability occurs at a thresh- 
old W* oc 1/a. If the dimensionless frequency detuning 
A| is less than a critical value A* ~ l/Q, the instability 
has a soft character, in that the amplitudes of the low fre- 
quency waves tend to zero at the threshold bifurcation 
line. Outside this range, the low frequency waves are 
characterised by the hard onset observed experimentally. 
Measurements (squares) are compared with theory (full 
line) in the main part of Fig. [5l The hard onset is accom- 
panied by a finite jump in subharmonic amplitude. These 
two regimes of behaviour are separated by critical points 
on the bifurcation line. The generation of subharmonics 
in the nonlinear oscillatory system found numerically ap- 
pears to be similar to the bifurcation of an anharmonic 
oscillator 2l|] . We estimated numerically that the crit- 
ical detuning parameter A* at T = 2.08 K is close to 
1CP 4 , and that the critical AC heat flux density W* is 
equal to a few mW/cm 2 , in good agreement with our 
observations. The experimental results shown in Fig. [T] 
can be understood as corresponding to the working point 
moving horizontally on the bifurcation diagram (Fig. [5l 
inset) into the (yellow) shaded region from a position just 
outside it, through the hard (blue) instability line. Al- 
though nonlocal (in w-space) interactions between waves 
at ujd and low frequency waves probably contributes to 
the subharmonics, due to the finite width of the low- 
frequency spectral domain, our numerical estimates show 
that the interactions of high-frequency (u ^> ujd) waves 
with subharmonics are relatively weak, consistent with 
our inference of two cascades. 

Our experiments have revealed an inverse wave energy 
cascade in a turbulent acoustic system. It is responsible 
for a substantial increase in wave amplitude correspond- 
ing to the formation of a set of huge low-frequency sub- 
harmonics. This instability develops through formation 
of isolated low-frequency waves of amplitude higher than 
that typical of waves around them. The latter can be 
considered as the acoustic analogue of the giant "rogue" 
waves that occasionally appear on the ocean and endan- 
ger shipping. Their origin lies in the development of a 
decay instability of the periodic wave, i.e. a similar mech- 



anism to that proposed [12j, |22| (modulation instability) 



to account for the creation of oceanic rogue waves [23j. 
Note that this mechanism differs from an alternative ex- 
planation proposed recently [lij that involves scattering 
of nonlinear waves on a continuous noisy background. 
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FIG. 5: (Color online) Dependence of the AC heat flux den- 
sity W at which the instability develops on the dimensionless 
frequency detuning A = (ujd — ui n )/uj n of the driving force fre- 
quency ujd from a cavity resonance u n . Numerical calculations 
(line) are compared with measurements (points) for driving 
at the 96th resonance. Horizontal bars mark the widths of the 
hysteretic region where second sound exists in a metastable 
state. Inset: bifurcation diagram showing regions of stabil- 
ity (unshaded) and and regions of instability (yellow shaded) 
against the generation of subharmonics. The soft instabil- 
ity occurs over the (orange) line between the (green) critical 
points at ±A*; outside them lies the hard instability; W* is 
the threshold value of the instability. 
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